Charge occupancy of two interacting electrons on artificial molecules — exact results 



A. Aharony 

School of Physics and Astronomy, Raymond and Beverly Sackler Faculty of Exact Sciences, 
Tel Aviv University, Tel Aviv 69978, Israel 

O. Entin-WohlmanS 

Centre for Advanced Studies, The Norwegian Academy, Oslo 0271, Norway 

Y. Imry and Y. Levinson 
Condensed Matter Physics, The Weizmann Institute of Science, Rehovot 76100, Israel 

(February 1, 2008) 

We present exact solutions for two interacting electrons on an artificial atom and on an artificial 
molecule made by one and two (single level) quantum dots connected by ideal leads. Specifically, 
we calculate the accumulated charge on the dots as function of the gate voltage eo, for various 
strengths of the electron-electron interaction U and of the hybridization between the dots and the 
(one-dimensional) leads 7. For 7 < f and 2(7 — 1) = eoo < eo < 0, there are no bound states. As 
eo decreases beyond eoo, the accumulated charge P in the two-electron ground state increases in 
gradual steps from to f and then to 2. The values P ~ represent an "insulating" state, where 
both electrons are bound to shallow states on the impurities. The value of P m 1 corresponds to 
a "metal", with one electron localized on the dots and the other extended on the leads. The value 
of 2 corresponds to another "insulator", with both electrons strongly localized. The width of the 
"metallic" regime diverges with U for the single dot, but remains very narrow for the double dot. 
These results are contrasted with the simple Coulomb blockade picture. 
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I. INTRODUCTION 

Much of the early work on quantum dots concentrated 
on large dots, which may contain many electrons, in the 
limit of very weak coupling between the dots and the 
leads, where one can employ the simple Coulomb block- 
ade picture. [Q In that limit, the energy cost of adding 
a single electron to the dot is of order of the excess en- 
ergy for the added charge. For a total of two electrons, 
this is of the order of the Coulomb repulsion between 
the two electrons. One expects that for a good contact 
between the dots and the leads this energy cost will be re- 
duced. However, the study of a general coupling strength 
presents a great challenge. The case of an 'open' dot, con- 
nected by quantum point contacts to a two-dimensional 
electron gas, Q has been analyzed using bozonization 
techniques and mapping of the Hamiltonian onto the two- 
and four-channel Kondo problem. || 

It has only recently become possible to also study small 
quantum dots, which have a small number of states and 
contain a small number of electrons. |^,^| Such a small 
dot, connected to external leads, is similar to a donor 
in a doped semiconductor: both may be modeled as an 
'impurity' connected to external leads. Q A set of such 
quantum dots, or an artificial molecule, can then be mod- 
eled by a set of such 'impurities'. In what follows, we 
sometimes interchange the terms 'dot' and 'impurity'. It 
is usually assumed that the electrons interact with each 
other only when they are on the same quantum dot, and 
behave as free electrons when they are on the leads. In 



what follows we therefore assume a contact interaction, 
which exists only on the dots, and present exact results 
for the case of two electrons. Given the difficulties in solv- 
ing the general problem, such analytical results (even for 
the most simple configruations) are helpful. They are 
particularly useful for nanostructures, where one might 
design controlled experiments. |^,|| 

We have recently reported on several exact results for 
two interacting electrons on a general number of dots 
A/", which are modeled as 'impurities' which have single 
electronic states: we presented a general scheme for find- 
ing the eigenenergies, and presented some results for the 
spectra of a single dot. Q We have also discussed the 
exact two-electron current through a single dot. || Here 
we generalize these results, with emphasis on the charge 
accumulated on each quantum dot and on its relation- 
ship with the Coulomb blockade picture. We then devote 
most of this paper to discuss the more complex case of a 
double dot. Such a double dot, with one state per dot, 
has recently been proposed as a possible candidate for 
the two-qubit entanglement required for quantum com- 
putation. H The case of two coupled quantum dots is 
also amenable to experiments. [[l0| 

In our earlier work 0, we showed that the spectrum 
and the wave functions of the two interacting electrons 
can be obtained in terms of the energy spectrum and the 
wave functions of the single-electron Hamiltonian. We 
reproduce these results in section II in a slightly differ- 
ent method, and use them to obtain new results for the 
charge occupancies on the quantum dots. The next two 
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sections are devoted to the study of specific configura- 
tions: a single dot, and a system made up of two dots, 
separated by a distance R. The single-electron spectra 
of these two configurations, required for the study of the 
two electron one, are discussed in the Appendix. 



II. TWO INTERACTING ELECTRONS - 
GENERAL SCHEME 

As has been demonstrated in Ref. M, the knowledge 
of the spectrum of the single-electron Hamiltonian is suf- 
ficient for deducing the spectrum and the wave functions 
of two interacting electrons, for any number A/" of dots. 
Basically, we start with the Hamiltonian 



H = H s 



H 



hit • 



(i) 



The spin-independent single-electron part TL ac involves 
site energies ti on the dots (i = 1, 2, A/") and zero on 
the lead sites, and also nearest neighbor hopping matrix 
elements t n ^ m which assume special values near the dots. 
This part is diagonalized by the eignenergies {e a } and 
the corresponding eigenfunctions {4> a (n)}. 

For simplicity, we assume that the two electrons in- 
teract only when they are both on the same dot i, with 
interaction energy C7(i) (though the method of solution 
can be extended for other types of interactions): 



(2) 



Using the single-electron eigenstates, the two-electron 
Hamiltonian takes the form 

^ = E £ a C L c aa + ^ E Uacbd W^T^T C cl Cd i ' 
acr i abed 

U acbd (i) = tf(i)#(i)#(i)&(i)6,(i). (3) 

Here, c\ a = 't'ain)^^ creates an electron in the state 
a with spin a. 

For such a contact electron-electron interaction, one is 
interested only in the singlet state of the two electrons 
(the energies of the two electrons in the triplet state are 
simply given by the non-interacting sums e a + e^). Wc 
hence write for the two-electron singlet wave function 



l*) = £x ab (£) C t T 4jo>, 



(4) 



ah 



where |0) is the vacuum and X ab = X ba . The Schrodinger 
equation 



H\V) = E\^) 



(5) 



then yields 



^(E-e a -e b )x db {E)cl,cl l \0) 

ab 

= E E x ab (E)u a , b , ab (i)cl^ bll \o). 



(0) 
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Multiplying this equation from the left by (0|o,"^c a "f 
gives 



x ab ( E )=j2i: 



U aba ' b ' i))Xa' b >{E) 

E - e a - e b 



(7) 



The simple form of the matrix elements of the contact 
interaction [see Eq. (J|)] allows us to rewrite Eq. (0) as 
a set of A/" linear equations: Defining the quantities 



Ai(E) =J2<p a (i)ct> b (i)X ab (E) 



(8) 



(which represent the amplitudes of for the singlet 
state with both electrons on site i, denoted by |i,i}), one 
arrives at 



A i (E) = ^U{i')G E (i,i ] i',i')A i ,(E), 



(9) 



in which Ge(1, i; i', i') is the two-particle Green's function 
of two non-interacting electrons, 



G E (n 1 ,n 2 ;n' 1 ,n' 2 ) = 



E - e a - e b 



(10) 



calculated at the impurity locations. |TT[ The determi- 
nant of Eqs. (||) gives the eigenenergies {E} of the two 
interacting electrons, and in particular determines the 
ground state energy, Eq. The coefficients X(E) are then 
obtained from Eq. (^), which can be rewritten as 



X ab (E) = J2 



U®ti (i)#(i)Ai(i^ 
E — ea — e& 



(11) 



Substituting this result into Eq. (|J), it is easy to check 
that 

I*) =J2 U ®ME) G E (n u n 2 ;i,i)\ ni ,n 2 }, (12) 



ni, n 2 ) is a singlet state with one electron at site n\ and 
the other at site n 2 . Note that Eq. (|j]) determines the 
Ai's only up to a multiplicative constant. This constant 
should be determined by the normalization of 1^}, i. e. 
from the condition ^2 ab \X ab (E)\ 2 = 1. Such solutions 
will be discussed in some detail below. 

The electronic states on a quantum dot are commonly 
probed by varying the gate voltages on the dots, repre- 
sented here by the ej's, and measuring the conductance. 
Alternatively, one may probe the total charge on the dots, 
by measuring the capacitance when the gate voltage is 
changed. Q The total charge on the dots, in the state 
is given by (in units of the electron charge, e) 
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^ = EE<*i c U^> 



AT 



(*lEE^(i)0 b (i)E c - c ^l*)- 



(13) 



i— 1 ab 

we obtain 



Using Eq. 

P = 2 E E ^(i)^(i)*a*c(£)*6c(£). (W) 



i afcc 

Alternatively, we note that 



P 



V — 



(15) 



This follows from first-order perturbation theory: writing 
E = (^\H\^f), the derivative with respect to t\ becomes 
(^/Ic^Cio-l^l/). In the following sections, we will use these 
results to discuss the ground state of simple quantum dot 
systems occupied by two electrons. 

An important technical advantage of the representa- 
tion of the two-electron spectrum in terms of the two- 
particle Green's function of the non-interacting system, 
is that the latter can be expressed in terms of the single- 
particle Green's function. The spectral representation of 
the single-particle Green's function, (^(n; n'), is 



g w (n;n') = ^ 



(16) 



u + i(-e a 

where ( — > + . In conjunction with Eq. ([lO]), one finds 

G E {ni,n 2 ; n[,n' 2 ) 
1 



2n 



dug E -uj(ni; ni)% w (n 2 ; n' 2 ) 

O 

dugs-u, {n x ; n[)g u (n 2 ; n' 2 ) , 



(17) 



where the last equality follows from the Kramers-Kronig 
relations. The relation ( |l7| ) is very useful in the detailed 
calculations of the ground state properties, because the 
single-particle Green's functions are relatively easy to 
find. We give in Appendix A the details of the single- 
particle Green's functions required for the dot configura- 
tions investigated in this paper. 



III. A SINGLE DOT ON A ONE-DIMENSIONAL 
WIRE 

We model a single dot on a one-dimensional wire by a 
single impurity, located at site 0, which has the on-site 
energy eo, and is coupled to two ideal one-dimensional 
leads, with the amplitude to for tunneling between the 



impurity and its nearest neighbors on the leads. Q Both 
these parameters are experimentally accesssible: eo mod- 
els the plunger gate voltage on the dot, while to is related 
to the transmittance (or conductance) of the barriers be- 
tween each dot and the leads, which can be varied by 
changing the gate voltages on these barriers. The cor- 
responding amplitudes between sites inside the leads are 
equal to —t. Although some aspects of the solution of 
this problem were discussed in Ref. 0, we present here 
an alternative derivation, which is more adapted to the 
calculation of the occupation P and to the more compli- 
cated case of two dots. 

As shown in Appendix A, (?|] the single-particle Hamil- 
tonian of such a system has none, one or two bound 
states, depending on the values of the "gate voltage" eo 
and the "hybridization" ig = 7 (energies are measured 
in units of t). For the sake of concreteness, we concen- 
trate on the region where there is only one bound state 
below the band, of energy ep < —2. This occurs for 
e < e o =2(7—1) (see Appendix A). 

In the simple case of a single 'impurity', Eq. (JoJ) re- 
duces to a single equation, and the eigenenergies {E} of 
the two interacting electrons are given by the solutions 
of 



- = G £ (0,0;0,0) 



E 

ab 



l0q(O)| 2 |^(O)| 2 
E - e a - e b 



(18) 



It is easy to deduce the beahvior of Ge as function of the 
two-electron energy E from this equation. Ge is negative 
for E < 2e/3, decreasing from to -co as E increases from 
-co towards 2ep. As E crosses this value, it jumps to 
+00 and then decreases. The value E = —2 + ep marks 
the beginning of the two-particle continuous band states: 
one electron is bound and the other is in the continuum. 

As discussed in Ref. JtJ, Ge is finite at E = —2 + ep 
in the thermodynamic limit of infinite leads, due to the 
vanishing of the band state wavefunction <pk (with en- 
ergy ek = — 2cosfc) at the impurity site i = for k = 0. 
This value of G determines whether there is or there 
is not a bound state of the two interacting electrons: 



When l/U < G-2+ep, then Eq. ( |18|) has no solution 
for E < — 2 + ep, and there is no doubly occupied bound 
state below the band. One of the electrons is then in 
a band state. The behavior of G at E = —2 + ep, as 
function of e , is plotted in Fig. |[ Since GE=—2+tn nas 
a maximum, G max , there is always a doubly occupied 
bound state (or an "insulator") for U < 1/G max . For 
larger U, the equation G_2+ £/3 = l/U has two solutions, 
eo,- and eo,+ (which depend on U). For eo,- < eo < eo,+ 
one has no doubly occupied bound state, and the ground 
state of the two electrons lies in the continuum (i. e. 
represents a "metal" ) . Then, as the on-site energy eo be- 
comes more attractive, the bound state of the two elec- 
trons re-appears. As seen from Fig. |l], e ,- diverges to 
-co when U — ► 00, and this re-entrance then disappears. 
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The region between eo,- and eo,+ becomes narrower as 
7 increases, and for finite U it always disappears above 
some critical hybridization "f c (which diverges to oo as 
tf-oo). @ 

These "insulator to metal" [|l2| transitions of the two- 
electron ground state, from being bound to being in the 
continuum and back, are reflected in the occupancy P of 
the dot in the ground state [see Eq. (|l~5|)]. To find P, 
we re- write the equation for the ground energy Eq in the 
form 



1 

U 



du{D 



E G -ui) 



\ — lfis 



(19) 



where we have used Eq. (|l7|), and the results for rel- 
evant for this geometry [Eq. (A7)]. The imaginary part 



appearing in this expression has a delta-function contri- 
bution coming from the single-electron bound energy, and 
the contribution arising from the band states. Separating 
these two, we have 



1 

U 



D E . 



G-E/3 

2 



- / duj{D 

K J -2 



Eg—u 



(20) 



where r(u>) = (dD u / dbj)^ 1 is the residue at the bound 
energy pole u — ep. It is now straightforward to differ- 
entiate all terms in ( p0| ) with respect to eo, and obtain 
P = dE G /de . 

We have solved Eq. (|2^) for E G , by calculating the 
integral numerically. We have then computed the varia- 
tion of the occupancy P as function of the gate voltage eo, 
and the results are shown in Fig. ^|, for a comparatively 
large hybridization, and in Fig. [|, for a small hybridiza- 
tion 7. Generally, P starts at at eo = eoo = 2(7 — 1), 
when the single electron bound state just moves below 
the band (with the inverse localization length Kp = 0, 
i. e. zero weight on the impurity). P then grows as eo 
decreases. As seen from Fig. ||, e ,+ is quite close to 
e o = 2(7 — 1), due to the steepness of G [G diverges 
to — 00 as eo | 2(7 — 1)]. Therefore, the first transition 
into the "metallic" phase, at eo,+ , occurs when the local- 
ization lengths of the two bound electrons are still quite 
large, and their weights on the impurity (and thus also P) 
are relatively small. The calculation of P in this regime is 
not easy, due to numerical problems related to the above 
mentioned steepness. In any case, P reaches values close 
to 1 somewhere inside the "metallic" phase, i. e. for 
£ o,- < e < e ,+. We note that in the first "insulating" 
phase, which appears at e ,+ < e < 2(7 — 1), both elec- 
trons are bound on very shallow states, hence the small 
value of P. Thus, it is not enough to know P in order 
to determine the transport nature of the system. As eo 
crosses below eo,_, into the second "insulating" phase, 
P gradually increases towards 2, reflecting the strongly 



bound state of the two electrons. This gradual increase 
becomes steeper as the hybridization 7 becomes smaller, 
and the width of the "metallic" single electron occupancy 
regime (of order eo.+ — eo,-) increases with increasing U . 
Both of these facts are in qualitative accordance with the 
Coulomb blockade picture (where usually the derivative 
of P with respect to the gate voltage has peaks whose 
width increases with the hybridization and whose inter- 
peak distance increases with U). [Qj|] In fact, the distance 
between the JV'th and the (N — l)'th peaks is usually in- 
terpreted as the energy cost of adding the TV'th electron. 
However, it is usually very difficult to obtain quantitative 
estimates for these quantities in that picture. Further- 
more, the similarity of our results to the simple Coulomb 
blockade picture is completely lost as 7 increases towards 
and beyond 7 C : the width of the "metallic" regime then 
shrinks, and the there is a continuous gradual increase of 
P from to 2. 

Returning to Eq. (|lj), we now observe that for a single 
impurity, the two-electron state is given by 



I*) = UA G E (n 1 ,n 2 ;0,0)\n 1 ,n 2 }, 



(21) 



where A is found from the normalization. One can now 
use Eq. (|j]) and the single electron Green's functions 
g UJ (n;0) to obtain For the bound ground state, the 
results show an exponential decay of the amplitudes as 
cither electron moves away from the impurity. 



IV. TWO DOTS ON A ONE-DIMENSIONAL 
WIRE 

Two dots are modeled by two 'impurities', connected 
to each other and to the outside by ideal linear leads. 
The presence of two impurities gives rise to up to four 
single-particle bound states. For simplicity, we consider 
two identical impurities, each having the same on-site 
energy eo, which are located at sites I and r, and are 
separated by a distance R (R > 2). Confining ourselves 
again to the configuration where the bound states appear 
only below the band, the first bound state appears when 
eo is smaller than eoo = 2(7—1), while the second appears 
only for R > R c , where 



R c = 27/(27- 2 -e ). 



(22) 



At fixed R, there exists a single bound state below the 
band only in the narrow regime 



2(7 - 1 - 7 /i?) < e < 2(7-1). 



(23) 



We also restrict ourselves to the regime with eo < 2(1—7), 
so that there are no bound states above the band (see 
Appendix A). 
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Assuming the on-site Coulomb interaction to be iden- 
tical on the two impurities, U(£) = U(r) = U, Eqs. (||) 
yield 

A e (E) = UG E (d)A e (E) + UG E (nd)A r (E), 

A r {E) = UG E (nd)A e (E) + UG E (d)A r (E), (24) 

where the labels d and nd stand for the diagonal and the 
nondiagonal elements of the matrix. By symmetry, 

G E (d) = G E (£, £;£,£) = G E (r,r;r,r), 
G E (nd) = G E (£, £■ r, r) = G E (r, r ■;£,£). (25) 

The eigenenergies of the two interacting electrons are 
given by the solutions of the two equations 



^=G E (d)±G E (nd) = G E , 

and the corresponding solutions obey 
Af(E) = ±At(E). 



(26) 



(27) 



It is instructive to rewrite these equations in terms of 
the single-electron wave functions, using Eq. ([lO]). In the 
symmetric molecule case, one can divide the solutions 
into even and odd single-electron wave functions, with 
(j> a {£) = ±</> a (7"). From Eq. (M) it now follows that 



G 



[MWbW±<i>a(r)<l>b(r)mWl(e) 



ab 



E 



(28) 



Thus, it is clear that contains only pairs of states 
where both a and b are even or odd, while G E contains 
only mixed combinations, where one state is even and 
the other is odd. It thus follows that the solutions of the 
two-electron problem divide into two separate families: 
the solutions of G^ = 1/U involve only even-even and 
odd-odd single electron states, while those of G E = 1/U 
involve only even-odd states: the coefficients in Eq. (^) 
will split into two separate families, associated with the 
different solutions of Eqs. (|2"6|). This can be easily seen 
by substituting Eq. @ into Eq. 

To discuss the two-electron energies, we need to ana- 
lyze the ^-dependence of G E . This depends on R: For 
R> R c , there exist two single-electron bound states, the 
even <j)p + and the odd <^g_. Thus, G^ decreases from oo 
to — oo as E increases from 2ep + towards 2e^_. There- 
fore, in this case the equation = 1/U always has a 
discrete solution, with E between 2ep + and 2e ( g_. In the 
same case, G E decreases from oo towards a finite value 
as E increases from ep + + ep- towards the bottom of the 
continuum —2 + ep+ (which contains both even and odd 
states). The new lowest even-odd state may thus be ei- 
ther "insulating" or "metallic" , depending on the sign of 
1/U — GZ 2 +c l}+ - However, the energy of this even-odd 
state is always above the lowest triplet energy, which is 



equal to the non- interacting value ep + + tp-. From our 
numerical calculations we observe that G^ is negative 
at this lowest triplet energy. Therefore, the lowest solu- 
tion of G^ = 1/U is the ground state of the two-electron 
problem, which is thus a singlet. In a way, this might 
have been expected: breaking the system into two parts, 
by removing a bond in the middle between the two dots, 
one ends up with separate "atomic" states on each side, 
each coupled to its own lead. Each side can then contain 
one electron with either spin up or spin down. However, 
switching on the hopping th between the two sides would 
lower the energy of the singlet state, similarly to the an- 
tiferromagnetic ground state of the Hubbard model; to 
lowest order in th, the "exchange" difference between the 
triplet and singlet states is of order t\/U. |l3|] It is in- 
teresting to note that in our case there exists a finite 
difference between the singlet and the triplet even in the 
limit U — > oo, since we find that G~^ +e is strictly neg- 
ative. It would be interesting to study generalizations of 
our simple model, e. g. including interdot Coulomb and 
exchange interactions, which would allow an interchange 
of the singlet and triplet ground states. |L4| 

The only chance to find a "metallic" ground state is 
thus for R < R c , when there exists only one single- 
electron bound state below the band. This limits the 
possible range of parameters to that in Eq. (^3|). In 
this regime, G E yields no doubly bound state, and G E 
yields one only if G^ 2+£/J+ < 1/U. Note that this regime 
becomes narrower (in terms of eo) as R increases. It 
is therefore interesting to find the borderline in the pa- 
rameter space, at which Gi 2 + C|3+ — 1/U. Inside this 
broderline, the ground energy of the two electrons is in 
the continuum, i. e. "metallic" . 

To calculate G E , we use Eq. ( p7| ) and the results ( A14) 
and (A15) of Appendix A: 



G E = G E (d) + Ge(nd) 



D 



E- 



1 



<j 

D+ .. D+ 



G E = G E (d) — G E (nd) 
1 




D 



E-u. 



(29) 



where are given by Eqs. ( |A15| ). 

We next separate the contributions of the bound ener- 
gies from the integrals, to find 

Gt = l(iM + G(R-R c )r~^ 




2 \ D 



D+ E-^ 



5 



— / dLu\ 



D 



— + ( 30 ) 



D , 



where r ± (cj) = {dD^j / 'duo) -1 are the residues at the 
poles. 

We have used Eq. ( |30| ) to solve the equation G± 2 + e)3 = 
0, which yields the borderline of the "metallic" regime in 
the limit U — * oo. The result for R — 2 is depicted by the 
dotted line in Fig. || The area enclosed inside this line 
represents the "metal" , where G~^ 2 + ei3+ > 0. For smaller 
U and for larger i?'s this area shrinks further. 

Figure [| highlights a major difference between the 
single-dot and the double-dot cases. In the former, the 
width of the "metallic" regime (in terms of eo) was equal 
to e ,+ — e ,_ , and at fixed 7 it increased with U, diverging 
to 00 for U — ► 00. Although this width was not equal to 
U, as assumed in the simple Coulomb blockade picture, 
it still resembled the qualitative features of that picture. 
In contrast, in the double-dot case this width is bounded 
by Eq. (p3|), and this bound is independent of U. There- 
fore, the width of the singly occupied "metallic" regime 
remains finite and small even when U — > 00. Basically, 
this happens because in the double-dot case, there exist 
two single-electron bound states. The level-repulsion be- 
tween these bound states then prevents the two-electron 
ground state from merging into the coninuum. In fact, 
we expect similar bounds on the Coulomb-blockade-like 
energy even for a single quantum dot, whenever the dot 
has more than a single bound state. 

As R increases, the two single-electron ground energies 
of the double dot become closer to each other and to 
the single-dot bound state energy. Since 2ep + < Eq < 
2ep-, it follows that the "insulating" ground energy of 
the two interacting electrons is almost independent of 
U. Moreover, as |eo| increases, the two single-particle 
bound energies approach one another, and practically we 
have only one, doubly-degenerate, single-particle bound 
energy. This behavior is shown in Fig. |^, for a very 
'open' dot. Similar effects arise when the hybridization 
is reduced: e/3+ and also become indistinguishable. It 
hence follows that independently of U, the ground state 
energy of the two interacting electrons is Eq ~ 2ep+ ~ 
2e / 3_. In such a situation, the charge accumulated on the 
quantum dot, P, will just follow the weight of the single- 
particle localized wave functions on the impurities. These 
have a 'smooth' behavior as function of eo (see Fig. ||). 
Consequently, the "Coulomb blockade" type behavior, 
which is obtained for the single- impurity dot (Figs. |^ 
and|^) is washed out. 

Finally, we comment on the two-electron wave func- 
tion in the ground state. At small R, when there exists 
only one single-electron bound state, this wave function 
is dominated by the even-even state in which both elec- 
trons occupy the state 4>p+. To leading order in 7 this 
term represents the Heitler-London molecular state, pro- 
portional to \t, I) + r, r) + \t, r) + |r, £), with equal weights 



to the electrons being on the same impurity or each elec- 
tron being on a different impurity. For large R, Eq is 
close to both 2ep+ and to 2e^_, and therefore the co- 
efficients X a b(EQ) will be dominated by the terms with 
a,b = /3+, j3+ or a,b = j3+, f3— , with roughly equal mag- 
nitudes for these two terms [see Eq. (0)]. Combining 
these two terms, it is easy to see that | \E r ) is then domi- 
nated by the atomic orbitals \£, r) + \r, £). Thus, our exact 
results interpolate nicely between these two leading ap- 
proximations, which are common in chemistry textbooks. 



V. CONCLUSIONS 

Our exact solution of the two-electron problem does 
not have the Fermi gases on the leads; it is limited to 
the 'canonical ensemble', with a fixed small number of 
electrons in the system. We do believe however that this 
exact solution for this simple case is still useful, in that 
it throws light on issues which sometimes remain unclear 
in approximate (and much more complicated) treatments 
of the real problem. 

Our solution may also be directly connected with e. 
g. the ionization of donors into the conduction band, 
as function of the system parameters. Simple models for 
this problem may involve N coupled one electron donors, 
or a single dielectronic donor. Our calculation shows that 
as the distance between a pair of single-electron donors 
decreases, then the on-site interaction U helps this ioniza- 
tion. This effect may well be an important ingredient for 
the real metal-insulator transition in some semiconduc- 
tors. Specifically, if each donor in the semiconductor has 
exactly one electron attached to it, then as the density 
of donors increases we expect pair of donors to combine 
into "molecules" which allow only one bound electron. 
The remaining electrons will move to the band, and the 
system will become "metallic" . 

The fact that even an onsite U can have effects which 
are so different from the naive Coulomb blockade model, 
should also be of interest. This is especially so for the 
double quatum dot case. As mentioned in the introduc- 
tion, the spectra of such systems are in principle address- 
able by transport and capacitance experiments. Usually, 
one looks at the contribution of the resonant states (lying 
in the continuum) of these dots to the conductivity, but 
the bound states will also contribute to the off-resonance 
transmission. The dependence of the average occupancy 
on the gate voltage is of interest both theoretically || 
and experimentally. 

Generalizations of this treatment to more realistic sit- 
uations, even for two electrons, are relatively easy to 
achieve. For example, the inclusion of an interdot inter- 
action V does not affect the need to solve only Af linear 
equations. As stated, such interactions may cause an in- 
terchange of the singlet and triplet ground states M , and 
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will certainly affect the magnetic exchange interactions 
between the electrons on different 'impurities'. 
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APPENDIX A: THE SINGLE-PARTICLE 
GREEN'S FUNCTION 

As discussed in the text, the two-particle Green's func- 
tion can be expressed in terms of the single-particle one. 
Here we derive the latter. 

The single-electron tight-binding type Hamiltonian, 
H sc , is given by 

n 

), (Al) 

n 

where we ignore the spin indices, since the single parti- 
cle Hamiltonian is spin- independent. Writing the single- 
particle Green's function, g, in the form 

g{n,n';t) = -*0(t)<[ct (t), (A2) 

we find, for the Fourier transform g u , the equation 

ujg^in, n') = S„^ + e n g w {n, ri) 

This equation is straightforwardly solved for the config- 
urations described in the text. 

1. The single impurity case 

For a single impurity with one-dimensional leads we 
have e n — 0, for n^O, and eo is the on-site energy of the 
impurity, t n>n ±\ = —t for n^O, ±1, and io,±i = £o- It is 
sufficient to consider Eq. (A3) for n' = 0. Then for any 
n ^ or ±1 that equation gives 



ug u (±n,0) = -tg u (±(n + l),Q)-tg u (±{n-l),Q). (A4) 

It is easy to convince oneself that, in the limit of infinite 
leads, the Green's function does not depend on the details 



of the boundary conditions. For n ^ one can therefore 
assume the solution 



g u (n,0) = C ul a^, 



and find from Eq. (A4) that 



J 



1/fltu 



(A5) 



(A6) 



Thus, for \ui/t\ < 2 we can denote a w = e lkuj , with uj/t = 

— 2cosfc w . For uj/t)< ± 2 we denote a u — ±e~ K ", and 
u>/t — =F2coshK w . 

The equation for n — ±1 now yields C w = 
— (to/t)g u (0, 0), and finally the equation for n = yields 



g u (0,0) = 1/D U , 

D tA) =u-e + 2(t 2 /t)a u 



(A7) 



From Eq. (|16|), the poles of g w (0, 0) give the eigenval- 
ues of the single electron problem, while the correspond- 
ing residues give the probability that an electron in a 
given state in on the impurity. The boundary between 
having or not having a bound state is easily found by 
setting a u = ±1 and uj/t = =p2 in the equation D u = 
(with the upper sign refering to a state below the band). 
Measuring energies in units of i, and denoting 



7 = (W*) 2 



(A8) 



we find that for eo < 2(7 — 1) there exists a bound state 
below the band (with < a u — e~ K P < 1), with a local- 
ization length I/k/j [which diverges to 00 at eo = 2(7 — 1)] 
and energy ep — — 2cos1i(k / 3). For e > 2(1 —7) there ap- 
pears a bound state above the band, with a localization 
length l/n a and energy e Q = 2cos1i(k q ). The two local- 
ization lengths are given by 



(A9) 



The weights of the localized wave function on the impu- 
rity (i.e., the residues of g u at the bound energies) are 
accordingly 



I^.MO)! 2 







de 



1 



27 



1 



(A10) 



A similar analysis gives the band of extended states, with 
energies efc = — 2cos/c. 

2. The two impurity case 

Here we consider a system with two impurities, which 
are separated by a distance R [R > 2). We denote the 
locations of the two impurities by £ and r, and assume 
t n ,n±i = ~t for n ^ £ or r, i„ )Ti ±i = t for n — £ or r, 
e„ = for n ^ i or r, and ee. r = eo. 
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Again, it is sufficient to consider g u (n,n') with n' = 
I, r. Referring to Eq. (|A3|), we assume a solution of the 
form 

g„(n,l)=C<(t)a£- e \, n<£, 
gUnJ) = C>(£)a^- r K n > r, 

g LJ (n,i)=A L0 {t)aZ + B UJ {l)aZ n , £ < n < r, (All) 



where a w is the solution of Eq. (A6). Writing Eq. ( |A3| ) 
for n = I — 1 and n = r + 1 gives 

C<(£) = -^g u (l,e) t C>(£) = - t j-g w (r,e). (A12) 

The other two coefficients, A^it) and B u {£) are found by 
using the equation for n — £ + 1 and n = r — 1. Then 



A u {t) = - 



to gu>{ZA)aJ ~ 9uj{r, £)a u 
t a^ R - aS 

t guj(r,£)ai -gu(ej)al 
t aZ R - a R 



(A13) 



Finally we write Eq. (A3) for n = £,r and obtain 



9u(r,r) = g u (£,£) 
9u,(£,r) = g u [r,£) 




(A14) 



where 



£o + — a u 



to a u> 



R+l 



,R-1 



± 



to a u 



eo + 7 ( a-u + 



a J 



t a u 
-l 



(A15) 



The single-particle bound energies are determined by 
the poles of the Green's functions, i.e., when van- 
ishes. Let us for simplicity confine ourselves to bound 
states below the band, with < a u < 1. Then D+ 
produces a bound state with energy e/3+ as long as 
eo < 2(7 — 1), as is the case for the single impurity 
configuration. However, the second bound state, e^-, 
coming from D~ , appears only at more negative eo, or 
(for fixed eo) when the distance between the impuri- 
ties, R, is large enough: eo < 2(7 — 1 — 7/i?) (solve 
D~ = with a u — * 1 and lu = —2). As R tends to 00, 
the two bound energies are approaching the same value, 
that of the bound energy of the single impurity system. 
We exemplify this behavior in Fig. |?], for 7 = 0.4 and 
eo = —1.5; the state with the higher energy appears only 
for R > R c ee 7/(7 - 1 - e /2) = 2.666.... The cal- 
culations presented in this paper are also restricted to 



eo < 2(1 — 7), so that there exist no bound states above 
the band. 

Generally, all the eigenstates of the problem divide 
into two subsets. Those which arise from D^j = obey 
gui(£, r)/g u {r, r) = ±1, and hence represent even (odd) 
solutions which obey 



gt{£-rn,t) = ±gt{r + m,i). 



(A16) 



In particular, one can associate these subsets with the 
even and odd single-electron wave functions, <f> a {£ — m) = 
±</>a(f + m)- The two bound states below the band, 
< / ) /3±("), thus correspond to the "bonding" and "anti- 
bonding" states of molecular chemistry. The antibonding 
energy e/3_ joins the band for inter-impurity distances be- 
low R r . 
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Raymond and Beverly Sackler Faculty of Exact Sciences, 
Tel Aviv University, Tel Aviv 69978, Israel. 
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FIG. 1. Ge(0,0;0,0) at E = -2 + e , as function of eo. 
The hybridization is fixed at 7 = 0.4. 
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FIG. 3. The occupancy of the single dot, as function of 
eo- The curves, from right to left, are for (7=2, (7 = 5, and 
= 10, respectively. The hybridization is fixed at 7 = 0.02. 
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FIG. 2. The occupancy of the single dot, as function of 
to- The curves, from right to left, are for (7=2, (7 = 5, and 
U — 10, respectively. The hybridization is fixed at 7 = 0.4. 



FIG. 4. eo _ 7 phase diagram for R = 2. The strai ght lines 
which meet at eo = —2 represent the bounds in Eq. (123), be- 
tween which there exists only one single-electron bound state 
below the band. Inside the dotted curved line, the system is 
"metallic", for U = 00. 
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FIG. 5. e/3+ and tp- as function of eo- Here R = 6, and FIG. 7. The two single-particle bound state energies below 

the dot is 'open': 7 = 0.98. the band, as function of the distance between the impurities. 
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FIG. 6. The weight of the single-particle wave functions in 
the bound states on the impurities. Here R = 6, and the dot 
is 'open': 7 = 0.98. 
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